Search for: All records

We report on the results of a simulation-based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, is relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative losses is included in the calculation. We systematically vary plasma beta (βm) in the flows, the strength of the cooling (Λ0), and the exponent (α) of temperature dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a ‘spine’ is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depend on initial flow parameters and provide a simple analytical model that captures the basic features of the flow.

 

We investigate the structure of the continental lithosphere by combining two approaches: a systematic survey of abrupt changes in seismic properties detected byP‐to‐Sconverted body waves and an integrated geophysical‐petrological inversion for temperature and density in the upper mantle. We refine the global thermo‐chemical model WINTERC‐G in eastern North America by including detailed regional information on the crust into petrological inversions and combine it with the upper mantle layering beneath eastern North America yielded by anisotropy‐aware receiver‐function analysis. Eastern North America's Archean, Proterozoic and Paleozoic lithospheres show an excellent agreement between the depth to the 1,300°C isotherm that bounds the lithosphere and the depth range where converted waves detect abrupt changes in seismic properties. Boundaries with these abrupt changes reside within the rigid mechanical lithosphere and are uncommon in the convecting mantle beneath it. The boundaries include both impedance increases and decreases with depth, as well as anisotropy changes, and must have developed over the course of the assembly and evolution of the lithosphere. In the asthenosphere below, such heterogeneities appear to have been largely mixed out by convection. The existence of abundant interfaces with diverse origin can account for the commonly observed scattered signals from within the continental lithosphere and presents an alternative to the end‐member concept of the mid‐lithospheric discontinuity as a ubiquitous feature with a uniform origin. Generally, we can define continental lithosphere as a region of conductive heat transport and steep geotherm that is characterized by pervasive internal layering of density, elastic moduli and texture.

 

Train-induced vibrations act as potential powerful high-frequency source for imaging subsurface with higher resolution than typical ambient noise interferometry. In this study, we present results of seismic interferometry applied on three days of railroad traffic data recorded by an array of seismographs along a railway in Dublin, Ireland. Our virtual shot gathers show significant surface and body wave energy that could be used for seismic interferometry. Reflection sections obtained with our interferometry approaches applied on selected time windows of train-induced vibrations is consistent with nearby borehole data and an active seismic profile. The consistency of the results given by these approaches confirms that train-generated vibrations represent a valuable source of signal for high-resolution subsurface imaging. Furthermore, our results show spurious arrivals that are due to the train geometry and also the cross-correlation approach that needs consideration for body wave interferometry studies.

 

We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive (“b-dot”) probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large ( R M > 10), the plasma flow advects a magnetic field proportional to the current at the load. This enables us to estimate the flow velocity as a function of time from the delay between the current at the load and the signal at the probe. The supersonic flow also generates a hydrodynamic bow shock around the probe, the structure of which depends on the upstream sonic Mach number. By imaging the shock around the probe with a Mach–Zehnder interferometer, we determine the upstream Mach number from the shock Mach angle, which we then use to determine the ion sound speed from the known upstream velocity. We use the sound speed to infer the value of [Formula: see text], where [Formula: see text] is the average ionization and T e is the electron temperature. We use this diagnostic to measure the time-resolved velocity and sound speed of a supersonic ( M S ∼ 8), super-Alfvénic ( M A ∼ 2) aluminum plasma generated during the ablation stage of an exploding wire array on the Magpie generator (1.4 MA, 250 ns). The velocity and [Formula: see text] measurements agree well with the optical Thompson scattering measurements reported in the literature and with 3D resistive magnetohydrodynamic simulations in GORGON. 

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

 

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period. 

Global variations in the propagation of fundamental-mode and overtone surface waves provide unique constraints on the low-frequency source properties and structure of the Earth’s upper mantle, transition zone and mid mantle. We construct a reference data set of multimode dispersion measurements by reconciling large and diverse catalogues of Love-wave (49.65 million) and Rayleigh-wave dispersion (177.66 million) from eight groups worldwide. The reference data set summarizes measurements of dispersion of fundamental-mode surface waves and up to six overtone branches from 44 871 earthquakes recorded on 12 222 globally distributed seismographic stations. Dispersion curves are specified at a set of reference periods between 25 and 250 s to determine propagation-phase anomalies with respect to a reference Earth model. Our procedures for reconciling data sets include: (1) controlling quality and salvaging missing metadata; (2) identifying discrepant measurements and reasons for discrepancies; (3) equalizing geographic coverage by constructing summary rays for travel-time observations and (4) constructing phase velocity maps at various wavelengths with combination of data types to evaluate inter-dataset consistency. We retrieved missing station and earthquake metadata in several legacy compilations and codified scalable formats to facilitate reproducibility, easy storage and fast input/output on high-performance-computing systems. Outliers can be attributed to cycle skipping, station polarity issues or overtone interference at specific epicentral distances. By assessing inter-dataset consistency across similar paths, we empirically quantified uncertainties in traveltime measurements. More than 95 per cent measurements of fundamental-mode dispersion are internally consistent, but agreement deteriorates for overtones especially branches 5 and 6. Systematic discrepancies between raw phase anomalies from various techniques can be attributed to discrepant theoretical approximations, reference Earth models and processing schemes. Phase-velocity variations yielded by the inversion of the summary data set are highly correlated (R ≥ 0.8) with those from the quality-controlled contributing data sets. Long-wavelength variations in fundamental-mode dispersion (50–100 s) are largely independent of the measurement technique with high correlations extending up to degree ∼25. Agreement degrades with increasing branch number and period; highly correlated structure is found only up to degree ∼10 at longer periods (T > 150 s) and up to degree ∼8 for overtones. Only 2ζ azimuthal variations in phase velocity of fundamental-mode Rayleigh waves were required by the reference data set; maps of 2ζ azimuthal variations are highly consistent between catalogues ( R = 0.6–0.8). Reference data with uncertainties are useful for improving existing measurement techniques, validating models of interior structure, calculating teleseismic data corrections in local or multiscale investigations and developing a 3-D reference Earth model.